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Decay of the cosmological constant by Hawking radiation as quantum tunneling 
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We calculate the emission rate of Hawking radiation from the cosmological horizon by quantum 
tunneling approaches. Using the Hamilton- Jacobi and the null geodesic methods, two typical ob- 
servations are obtained. First, the spectrum of radiation is not strictly pure thermal. Second, the 
value of the cosmological constant decreases. Our calculation is different from other similar works 
about Hawking radiation from the cosmological horizon in that the horizon does not shrink but 
expands by taking into account the back-reaction effect. We show that when a positive frequency 
particle is materialized near the cosmological horizon, the contribution of a negative frequency part- 
ner to the background geometry lowers the value of the cosmological constant. We also touch upon 
thermodynamics of de Sitter space and the quantum creation of universe. 

PACS numbers: 04.70.Dy, 04.60.-m, 97.60.Lf 



I. INTRODUCTION 

Understanding of black hole physics is believed to be 
a hint for the quantum theory of gravity. Among others, 
Hawking radiation from black holes is one of the most 
striking effects that is known to arise from the combina- 
tion of quantum mechanics and general relativity. There 
are several derivations of Hawking radiation. Hawking's 
original derivation, which calculates the Bogoliubov co- 
efficients between in and out states for a body collapsing 
to a black hole, is the most direct and physical [l], 
Derivations based on Euclidean quantum gravity are very 
straightforward Path integral derivation by Hartle 
and Hawking uses the Kruskal extension of spacetime and 
the physical time is analytically continued to the imag- 
inary time. Then Hawking temperature of the space- 
time is equal to the imaginary time period of a semiclas- 
sical propagator [4;]- These derivations, however, treat 
the background geometry fixed, and energy conservation 
is not enforced during the emission process. Evidently, 
a black hole loses its mass when particles are emitted. 
Hence, conventional approaches to derive Hawking ra- 
diation can not describe more realistic mass decreasing 
phenomena. 

Recently, Parikh and Wilczek derived Hawking radia- 
tion as a quantum tunneling process where energy con- 
servation plays a fundamental role Their tunneling 
picture is based on a pair creation of particles near the 
event horizon with zero total energy. If a negative energy 
particle is closer to the center of the black hole than a 
positive energy particle, the positive energy particle will 
see, by virtue of Birkhoff 's theorem, as the horizon move 
inward and may find itself on the outside of the horizon. 
Then, if the positive energy particle is radiated away from 
near the outside of the contracted horizon to infinity, the 
mass of the black hole decreases because the negative 
energy partner remains inside the horizon of black hole. 



As a consequence, the horizon contracts from its original 
radius to a new smaller radius. In this process, the clas- 
sical one particle action becomes complex and gives the 
WKB amplitude an imaginary part. The emission rate 
r of Parikh and Wilczek provides the rate with a leading 
order correction arising from the loss of black hole mass 
which corresponds to the energy carried by the radiated 
quantum. The result is written down simply as 



P oc exp (ASbh) ■ 



(1) 
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where ASbh is the change of Bekenstein-Hawking en- 
tropy of the black hole between before and after the par- 
ticle emission. For the low energy limit, eq.{T]) becomes 
r oc exp (—(3E), where fi is the inverse Hawking temper- 
ature and E is the energy of the emitted particle. This 
is the usual Boltzmann factor. These facts indicate that 
the spectrum of black hole radiation is not exactly pure 
thermal. We call this method to calculate the rate T 
the null geodesic method. There have been considerable 
efforts to generalize this work to those of various black 
holes and the similar results which support Parikh and 
Wilczek's tunneling picture have been derived M , H , IE H, 
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The other approach to the tunneling from a black hole 
is the Hamilton-Jacobi method considered by Angheben 
et al.[13|- I n this method, the action of a emitted parti- 
cle satisfies the Hamilton-Jacobi equation. Solving this 
equation, one finds that the action has an imaginary part. 
This imaginary part arises from a classically forbidden 
path for the emitted particle. Due to vacuum fluctua- 
tions near the horizon, a pair of particles can be created 
inside the horizon. If the positive energy particle tunnels 
out the horizon to infinity, it is observed as a Hawking ra- 
diation by a distant observer. In order to assure this pro- 
cess, it is required that a positive energy particle moves 
across the horizon from inside to outside. This path is 
obviously a classically forbidden path. This process, how- 
ever, can be interpreted as a tunneling process quantum 
mechanically. We call this method the Hamilton-Jacobi 
method. Several investigations along this work have also 
been done and led to the correct Hawking temperatures 
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The difference of above two methods lies on the part 
where the classically impossible process occurs. The null 
geodesic method treats particles classically and the back- 
ground geometry dynamically, namely, quantum mechan- 
ically. It is assumed that particles follow the geodesies 
and move classically. This is the reason we call it the 
null geodesic method. This method has the classical 
impossibility in the background geometry because the 
horizon of black hole can not contract classically. On 
the other hand, the Hamilton-Jacobi method treats the 
background geometry classically, namely, as a fixed back- 
ground, and the particle which tunnels out the horizon is 
treated quantum mechanically. In the Hamilton-Jacobi 
method, the particle does not follow the geodesic but 
moves the classically forbidden path, namely, across the 
horizon path. 

In this paper, we study Hawking radiation from the 
cosmological horizon by using the Hamilton-Jacobi and 
the null geodesic methods. This issue has been studied 



by using the former method in ref. 2J| . But their result 
seems to have the wrong sign. The Hawking temperature 
with respect to the cosmological horizon has been given 
by eq.(4.8) of ref.[24|]. In the vanishing black hole mass 
limit or empty de Sitter space limit, rc — > I, their inverse 
Hawking temperature /3c reduces to —2irl, which is obvi- 
ously negative value. Since negative temperature seems 
to be unphysical, probably something may be wrong. 
In ref.[2J|, the detailed calculation is not given. Thus, 
we shall reconsider the Hawking temperature of the cos- 
mological horizon by using the Hamilton-Jacobi method. 
This is the first motivation of our paper. 

The calculation of the emission rate from the cosmolog- 
ical horizon by the null geodesic method also has been 
studied in refs.[l8l. [l9| for empty de Sitter space. At 
first sight, one may think that their calculation is con- 
sistent and correct. Although they have considered the 
process in which the cosmological horizon shrinks by tak- 
ing into account the particle's self-gravitation, this con- 
tradicts the original idea of the null geodesic method. In 
this method, it is assumed that the positive energy parti- 
cle follows its geodesic and moves the classically allowed 
path. If the cosmological horizon shrinks, the particle 
must move a classically forbidden path in order to materi- 
alize as Hawking radiation. In this paper, we consider the 
different tunneling process from the cosmological horizon 
using the null geodesic method and derive the new result 
which is different from refs.[l8|, [n|. This is the second 
motivation of our paper. In the null geodesic method, 
we consider the tunneling process that the cosmologi- 
cal horizon expands and the positive frequency particle 
moves along the classically allowed path. Although our 
tunneling process is different from that of refs.[18|, [l9(, it 
matches with the original idea of Parikh and Wilczek's 
tunneling picture Q. Our result is consistent with the 
underlying unitary theory, that is, the emission rate is 
written down in the form of eq. ([T]) . Furthermore, our re- 
sult indicates that the value of the cosmological constant 



decreases due to energy conservation |3fj |. 

The organization of this paper is as follows. In sec- 
tion HH we calculate the emission rate from the cosmo- 
logical horizon by using the Hamilton-Jacobi method. 
Although our result appears to be the unfamiliar ex- 
pression, eq. (|15[) . it is interpreted as a usual Boltzmann 
factor. Section IIIII is devoted to the explanation of 
the difference between our tunneling picture and that 
of refs.[l8[ [l9j . The emission rate by the null geodesic 
method is calculated in section HVl Finally, in section fVl 
we discuss our result of this paper. 

Throughout this paper, the metric signature adopted 
is ( — , +, +, +). The use has been made of natural units, 
namely ?i = c = G = las well as k = 1 . 



II. HAMILTON-JACOBI METHOD 

In this section, we calculate the emission rate from 
the cosmological horizon by using the Hamilton-Jacobi 
method. In order to explain its derivation in detail, we 
compare the calculation of the emission rate from the cos- 
mological horizon with that of Schwarzschild black hole. 
It is shown that the emission rate from the cosmological 
horizon does not have a familiar form but it is interpreted 
usual Boltzmann factor. 

We consider a static spherically symmetric spacetime 
described by the metric 



ds z 



-f(r)dt 2 + r\r)dr 2 + r 2 dil 2 , 



(2) 



where dfl 2 is the metric of two dimensional unit sphere. 
This metric includes the Schwarzschild solution and the 
de Sitter solution as special cases. In this background, a 
massless scalar particle follows the Klein-Gordon equa- 
tion V M V''0 = 0, where V M represents the covariant 
derivative with respect to the metric ([2]). We look for 
the solution having the form <f> cx exp (il + ■■■), where / 
denotes the classical action of the particle and dots repre- 
sents the higher order terms in h. If this classical action 
has an imaginary part, the emission rate T is given by 



r oc exp (—2 Im /) . 



(3) 



In the following, we calculate the emission rate along 
the similar line on ref.[3l| and then we check our result 
by using the proper distance approach which has been 
considered originally in ref.[13|. Note that in ref.[3l| the 

scalar particle is defined as cx exp (— il -\ ) so that 

some signs arc different from ours. 

Inserting 4> oc exp (il + ■ • • ) into the Klein-Gordon 
equation, one finds the following equation to leading or- 
der in h, 



g^d^ldj = 0. 



(4) 



This is the relativistic Hamilton-Jacobi equation for the 
classical action of a massless particle in curved spacetime. 
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By the symmetry of the background spacetime, we can 
use the separation of variables as follows, 



I = -ut + W(r), 



(5) 



where we have ignored the angular part because we con- 
sider the s-wave only in this paper, to is a constant which 
corresponds to the frequency of the particle. From eqs. ((4j) 
and |5|), the equation for W(r) becomes 



o 



and the solution of this equation is expressed as 

dr 



W{r) = ±uj 



fir) 



(6) 



(7) 



where upper (lower) sign represents that the particle is 
the outgoing (incoming) wave. 

Now we consider the case of Schwarzschild black hole, 
in which the function f(r) in the metric ([2]) is written as 
f(r) = 1 — rii/r. Of course, t\h = 2M is the radius of 
the horizon. For the Schwarzschild case, eq.® becomes 



W(r) = ±uj [r + r H 



dr 



(8) 



There is a pole at r = Th hi the r integral. In order to 
avoid the pole we use Feynman's prescription and replace 
r — m with r — m — *£■ This yields 



W{r) = ±uj 



iirrH + rjj P 



1 



r H 



dr 



(9) 

where a is an integration constant and P() denotes the 
principal value. We are interested in the imaginary part 
of the classical action. From eqs.® and ©, they are 
given by 



Im / = ±nujrH + Im a. 



(10) 



Incoming particles can fall behind the horizon along clas- 
sically permitted trajectories, with a capture cross sec- 
tion of order of the horizon area. Hence, the classical ac- 
tion for incoming particles must be real. This determines 
the imaginary part of a constant a, Im a = ttluth- There- 
fore the emission rate from Schwarzschild black hole is 
written as 



T oc exp (SttMui) . 



(11) 



This is nothing but the Boltzmann factor with the tem- 
perature Th — 1/8ttM which is the Hawking tempera- 
ture of the Schwarzschild black hole. 

Note that in order to obtain the emission rate we have 
used the boundary condition for the incoming particle. 
This is equal to considering the ratio between the particle 
emission and absorption. It is known that semiclassically 
this ratio is given by P(emission) = e~^ B P(absorption), 



where (3 is the inverse Hawking temperature and E is 
energy of the emitted particle (ij. 

Along the similar line, we consider the de Sitter case, 
namely, the emission rate from the cosmological horizon. 
For de Sitter spacetime, the function f(r) is given by 
f(r) — 1 — r 2 /r 2 H , where tb — I = \/3/A is the radius 
of the cosmological horizon and A is the cosmological 
constant. In this case, eq.([7|) has the form 



W{r)= T u J r f 



\og(r + r H ) + J 



dr 



a. (12) 



In order to avoid the pole, we again use Feynman's pre- 
scription and replace r — ry with r — r# + it instead 
of r — rn — it- This is because the positive frequency 
particle is the incoming wave for the radiation from the 
cosmological horizon. Although the sign of it is reversed 
in comparison with that of the Schwarzschild case, the 
positive frequency particle decays in time for both the 
Schwarzschild and de Sitter cases by these choices. Then, 
ea.(fl2]) yields 



W(r) = tw 



2 
P 



— in — log (r + ru) 



1 



dr 



(13) 



a, 



J-r H/ 

and the imaginary part of the classical action becomes 
Jml = ±— ujrn + Ima. (14) 
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For the de Sitter case, since outgoing particles can move 
away beyond the horizon along classically allowed paths, 
their actions are real. This gives Ima = —nujrH/2. 
Therefore the emission rate from the cosmological hori- 
zon becomes 



r tx exp (2ttIlo) . 



(15) 



If one identifies 0c — —^ttI as the inverse Hawking tem- 
perature of the cosmological horizon, the result of ref. [24[ , 
r oc exp (— 0coj), is reproduced. Since negative tempera- 
ture seems to be unphysical, however, we propose a differ- 
ent interpretation of the emission rate (TT5]) . Before doing 
it, we shall derive the emission rate (|15[) by an another 
approach in which the proper spatial distance is used. 
As Angheben et al. have used originally the proper spa- 
tial distance to calculate the emission rate |24| , we check 
whether the replacement r — rn — > r— rjj+ie is consistent 
with the proper distance approach. 

For the metric the proper distance a is defined by 



dr 



Using this proper distance, eq.([7]) is written as 



(16) 



(17) 
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Since, near the horizon, function f(r) is expanded as 
f(r) ~ f'(r H )(r — r H ) + • • • , the proper distance becomes 



-.y/r - r H . 



Then W(r) has the following form, 

2lu f da 



W{r) = ± 



\f'{r H )\ 



(18) 



(19) 



First, we shall consider the Schwarzschild case where 
|/'(rff)| = 1/2M. To avoid the pole at a = 0, we use 
Feynman's prescription a — > a — ie. Then the imaginary 
part of the classical action is given by 



Im I = ±2ttM w + Ima, 



(20) 



where a is a constant. Here we have used the fact that 
the integral contour is transformed from a semicircle to 
a quarter circle by coordinate transformation r — > a (32| 
which is given by eq. ([T5|) . This contour yields the factor 
in/2 on the imaginary part rather than in. Equation 
(|20f is the same with eq.tfTOjl so that the emission rate 
from Schwarzschild black hole is given by eq. (jTTJ) . This is 
consistent with the previous result for the emission rate. 

Next, we consider the de Sitter case. One may ex- 
pect that if the same calculation with the Schwarzschild 
case is done, the correct result will be derived. There is, 
however, the difference between the Schwarzschild case 
and the de Sitter case with respect to the definition of 
the proper spatial distance. For the Schwarzschild case, 
from cq. (|16jl the proper distance is given by 



a = yj r(r — 



rn) 



rn , 
T log 



+ y/r - r H 



if - y/r - r H 



(21) 



where we have chosen an integration constant as a = at 
r = rn- This proper distance is positive definite as long 
as r > rn- On the other hand, for the de Sitter case, if 
the same definition, eq.(fl6]). is used, one finds 



a = r H 



arcsm 



(22) 



As before, an integration constant has been chosen as 
a = at r = rjj. This proper distance is negative definite 
when r < rn- For the Schwarzschild case, the outgoing 
wave which has the positive frequency propagates toward 
the a increasing direction. On the other hand, for the de 
Sitter case, the incoming wave which has the positive 
frequency propagates toward the a decreasing direction. 

To avoid the sign ambiguity of the ie prescription, we 
use the same relation between the direction of motion of 
the positive frequency particle and the direction of a axis 
for both the Schwarzschild and de Sitter cases. In order 
to make a positive definite, we define the proper distance 



dr 



(23) 



for the de Sitter case instead of eq. ([To 
distance becomes 



a = r H 



7r i r 

arcsm — 

2 W 



Then the proper 



(24) 



This proper distance is positive definite when r < th 
and measures the proper distance from the cosmologi- 
cal horizon rn to a certain point r (r < rn)- By this 
definition, the incoming wave which has the positive fre- 
quency moves toward the a increasing direction as the 
Schwarzschild case. 

Note that the minus sign in eq. (|23|) arises from the fact 
that the upper bound and the lower bound of the integral 
is reversed. For the Schwarzschild case, in eci. tfil))) . it is 
implicitly assumed that the upper bound is larger than 
the lower bound. On the other hand, for the de Sitter 
case, we measure the proper distance from the horizon rn 
to a certain point r (r < rn ) so that the upper bound is 
smaller than the lower bound. To keep the larger radius 
as the upper bound of the integral, we have inserted the 
minus sign in eq. (|23[) . 

One can obtain the imaginary part of the classical ac- 
tion from eq. (|19[) with the replacement a — * — cr, but 
the expression for W(r) remains eq. (|19p . where the up- 
per bound is larger than the lower bound of the integral. 
Next, we use Feynman's prescription in order to avoid the 
pole. Since we have expressed the action by the proper 
distance, the positive frequency particle moves toward a 
increasing direction for both the Schwarzschild and de 
Sitter cases. So we need not now take into consideration 
the direction of motion when Feynman's prescription is 
used. In other words, for both spacetimes, we can use 
the replacement a — > a — ie equally. Consequently, the 
imaginary part of the action becomes 



Im7 = ±—ujrH 
2 



Im a, 



(25) 



where a is a constant. This is the same expression with 
eq. (fT4")) so that the emission rate from the cosmological 
horizon is given by eq. (|15p again. 

Finally, we propose a new interpretation of the emis- 
sion rate (fT5|) . As mentioned above, if one identifies the 
inverse Hawking temperature of the cosmological horizon 
as (3c = —2itI, the emission rate is written as 



T tx exp (-ficoj) , 



(26) 



which appears to be a usual Boltzmann factor and this 
agrees with the result of ref.[24|. Negative temperature, 
however, seems to be unphysical for the emission rate. 
Rather, it is natural to consider that temperature is al- 
ways positive. In this paper we stand this point of view 
and identify flc = 2itl as the inverse Hawking temper- 
ature. Then the Hawking temperature of the cosmolog- 
ical horizon is, of course, given by a well known result 
Th = l/2irl [331 ] . Using this temperature, the emission 
rate (fT5|) becomes T cx exp (/3coj). We write it as follows, 



T tx exp (-P C E) 



E 



(27) 
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In order to explain the physical meaning of eq. (|27[) , we 
shall recall the Schwarzschild-de Sitter solution which is 
represented by the metric (J2|) with 

, \ 2m r 2 . , 

m = i- — -p. (28) 

Hereafter we distinguish the terms "energy" and "mass" . 
The Schwarzschild-de Sitter solution approaches the 
usual Schwarzschild solution in the limit I — > oo so that 
the parameter m represents the mass of black hole. On 
the other hand, it is known that the energy of this sys- 
tem is given by E — —m [34[ . That is, positive "mass" 
is measured as negative "energy" in de Sitter space. 

Similarly, we shall distinguish the terms "energy" and 
"frequency" when we discuss the problem of the particle 
emission from the cosmological horizon in empty de Sitter 
space. If a emitted particle has positive "frequency" uj, 
it will be measured as negative "energy" E = — uj. If 
this particle collapses to a black hole near r = 0, then 
the metric after the collapse will be described by eqs.((2|) 
and (|28p with m = uj. This spacetime has the black hole 
mass uj and the energy E = —lu. It is reasonable to 
consider that this collapse will occur naturally because 
the total energy is conserved before and after the collapse. 
From these facts, we can recognize that for a particle 
with the positive frequency uj, its energy is expressed as 
E = —uj. Then the emission rate is written by means of 
this energy as eq. (f2"T)) . As a result, we find that eq. (fT5|) 
or (f27|) represents the usual Boltzmann factor with the 
inverse Hawking temperature (3c = 



III. NULL GEODESIC METHOD 

In this section, we consider the particle emission by the 
null geodesic method. Here we stress the fundamental 
difference between our approach and others. The emis- 
sion rate is calculated concretely in the next section. 

To begin with, we review the original idea of quantum 
tunneling by Parikh and Wilczek [5[. Their tunneling 
picture is based on a pair creation of particles near the 
horizon with zero total energy. If a negative energy parti- 
cle is closer to the center of the black hole than a positive 
energy particle, the positive energy partner will see the 
horizon shrink and may find itself on the outside of the 
black hole. The contribution of the negative energy parti- 
cle to the background geometry makes the event horizon 
contract from the original radius to the new smaller one. 
If the positive energy particle lies on the outside of the 
new horizon and propagates classically to infinity, it is 
observed as Hawking radiation. 

To describe across-horizon phenomena, it is necessary 
to choose coordinates which are not singular at the hori- 
zon. A particular suitable choice is Painleve coordinates. 



n = 2M 




r 



FIG. 1: Tunneling picture from Schwarzschild black hole. The 
event horizon contracts from the initial radius n — 2M to 
the final radius r/ = 2{M — uj). A pair creation occurs at 
n — e from which the outgoing positive energy particle starts 
to move along its geodesic. The particle materializes as a 
real particle at r/ + e, from which it propagates classically to 
infinity. 

For Schwarzschild black hole, it is written as follows 

ds 2 = - (l - dt 2 + 2^-dtdr + dr 2 + r 2 dQ 2 . 

(29) 

Then the radial null geodesies are given by 

,-±1-^, (30) 

with the upper (lower) sign corresponding to the outgo- 
ing (incoming) geodesies, under the implicit assumption 
that t increases toward the future. These equations are 
modified when the particle's self-gravitation is taken into 
account. If the total energy is fixed by energy conserva- 
tion, we should use eqs. ([29|) and l[30|) with M — > M — uj, 
where uj is the energy of the emitted particle. 

The null geodesic method assumes that the positive 
energy particle follows the classical orbit, namely, its 
geodesic. When a pair creation occurs inside the hori- 
zon, one finds, from eq. (f30|) . both incoming and outgoing 
particles move toward r = 0. Hence, particles can not 
escape from the black hole. In the idea of Parikh and 
Wilczek, one treats the background geometry dynami- 
cally. As mentioned above, the contribution of the neg- 
ative energy particle makes the horizon contract. This 
is the classically impossible process. Classically, black 
hole can increase the area according to Hawking's area 
theorem, but can not decrease it or radius. In the null 
geodesic method, this impossibility is not valid due to 
quantum mechanical effect by which the particle emis- 
sion can occur. By energy conservation, if the particle 
which has the positive energy ui is emitted, the radius 
of black hole must contract. The authors of ref.[f| ex- 
press these phenomena in their paper as follows. "Al- 
though this radially inward motion appears at first sight 
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to be classically allowed, it is nevertheless a classically 
forbidden trajectory because the apparent horizon is it- 
self contracting" and "the outgoing particle starts from 
r = 2M — e, just inside the initial position of the hori- 
zon, and traverses the contracting horizon to materialize 
at r = 2(M — ui) + e, just outside the final position of the 
horizon" . This situation is drawn in FigfT] There are two 
key points in the null geodesic method. First, the positive 
energy particle moves along a classically allowed path or 
its geodesic. Second, the contracting horizon passes the 
positive energy particle from behind. Therefore the ra- 
dius of the initial horizon must be larger than that of the 
final horizon of the black hole. 

Along the similar line, we consider the tunneling from 
the cosmological horizon. The line element of de Sitter 
space in Painleve coordinates is expressed as follows [l8[ , 

ds 2 = - (l - ^ dt 2 - 2 T -dtdr + dr 2 + r 2 dfl. (31) 

The radial null geodesies are given by 

r=y±l, (32) 

where upper (lower) sign represents the outgoing (incom- 
ing) geodesies. When the particle is beyond the horizon 
(r > I), both outgoing and incoming geodesies corre- 
spond to increasing r and the particle never crosses the 
horizon classically. This is the important fact in order 
to apply the null geodesic method to the cosmological 
horizon. 

Suppose that the pair creation occurs outside the hori- 
zon. If the incoming particle starts from r = I + e, just 
outside the initial position of the horizon, it moves to- 
ward r > direction according to its geodesic. In order 
to materialize at the point just the final position of the 
horizon, the cosmological horizon must expand. So the 
phenomenon drawn in Figf2]is required to emit the par- 
ticle from the cosmological horizon. Hence the contribu- 
tion of the negative frequency particle to the background 
geometry must make the cosmological horizon expand. 

Meanwhile, the quantum tunneling from the cosmolog- 
ical horizon has been already studied by Parikh [18[ and 
Medved [HI for empty de Sitter space. They, however, 
have assumed that the background geometry changes 
from empty de Sitter space to Schwarzschild-de Sitter 
space. If this happens, the radius of the cosmological 
horizon contracts (r/ < r%). In this case, if the positive 
frequency particle follows its geodesic, the tunneling can 
not occur, because the horizon does not pass the particle 
from behind. Both authors of refs.[l8|, [lj|, further, have 
assumed that the pair creation occurs just outside the ini- 
tial horizon and the positive frequency particle tunnels 
out just inside the final horizon. Indeed, Parikh claims in 
ref.[l8| as follows. "We expect r\ to correspond roughly 
to the site of pair creation, which should be slightly out- 
side the horizon" and "We expect 77 to be a classical 
turning point, at which the semiclassical trajectory (i.e. 
instanton) can join onto a classical-allowed motion. This 



r f = I + cj 




r 



FIG. 2: Tunneling picture from the cosmological horizon. The 
horizon expands from the initial radius r% = I to the final 
radius 77 = I + u (see section ITVl for the details of the final 
radius r/). A pair creation occurs at ri + e from which the 
incoming positive frequency particle starts to move along its 
geodesic. The particle materializes as a real particle at r/ — e, 
from which it propagates classically inward the origin (r = 0) . 



r, 




FIG. 3: Tunneling picture which has been considered in 
refs.[l8j, 19]. The cosmological horizon contracts from the ini- 
tial radius n to the final radius r/, because the background 
changes from empty de Sitter space to Schwarzschild-de Sit- 
ter space. A pair creation occurs at n + e. The path from 
rj+etor/ — e is the classically forbidden one. Thus, it is not 
allowed to use this path in the null geodesic method. 



must be slightly within the horizon." Under the assump- 
tion that the cosmological horizon shrinks, this situation 
is drawn as Figl3] In this process, obviously, the emitted 
particle does not follow the classically allowed path be- 
cause all particles move toward r > direction for r > I. 

From these facts, we think that the calculation of the 
tunneling rate from the cosmological horizon studied in 
refs.JHI, [l!| may be wrong or does not match with the 
original idea of quantum tunneling constructed in ref. [j| . 
By these reasons, in the next section, we investigate the 
emission rate from the cosmological horizon using the null 
geodesic method based on FiglJJ The fundamental differ- 
ence between ours and refs.fla, [l9j lies on the change of 
the horizon radius and the path of the positive frequency 
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particle. 



IV. THE EMISSION RATE 



In this section, we calculate the emission rate from 
the cosmological horizon based on the consideration of 
previous section and FigO 

The imaginary part of the action of an s-wave incoming 
particle which crosses the horizon can be written as [j| 



Im / = Im 



Pr 



Im 



dp r dr. (33) 



Using Hamilton's equation, we change the variable from 
momentum to energy, then we have 



Im / = Im 



dH , 
— —dr. 

r 



(34) 



Here we must consider two problems. The first prob- 
lem is what is the total energy or Hamiltonian of de Sit- 
ter space. The other is how to include the back-reaction 
effect by particle's self-gravitation. We can estimate the 
variation of Hamiltonian dH easily. In section|TTl we have 
given the energy of the emitted positive frequency parti- 
cle as E — —uj. Since the pair creation occurs with zero 
total energy, the negative frequency partner should have 
the energy E = +uj. Thus, the variation of Hamiltonian 
is given by dH = +duj. Meanwhile, for the change of the 
horizon radius, we can also have the simple estimation 
as follows. First, temperature and entropy of de Sitter 
space is Th = l/2nl and S — nl 2 . Since de Sitter space 
is one parameter spacetime as Schwarzschild spacetime, 
one can deduce that temperature, entropy and energy of 
the system satisfies the relation 



dE = T H dS. 



(35) 



For de Sitter space, the right hand side of eq. (|35|) yields 
TndS = dl. Hence, the total energy of de Sitter space 
should have the form E = 1 + const., where constant term 
is not important for latter considerations. Thus we can 
use H = I as the Hamiltonian. When the contribution 
of the negative frequency particle to the background ge- 
ometry is taken into account, the energy +uj is added 
and one has W = I + uj as the modified Hamiltonian. 
Therefore, the effect of back-reaction makes the horizon 
expand from the original radius r; = I to the new larger 
one r/ = I + uj. This is the expected result, because the 
expansion of the radius of the cosmological horizon does 
not contradict with the tunneling picture of Fig[2] 

Although the fact that the contribution of the nega- 
tive frequency particle increases the energy of the system 
appears at first sight to be curious, this can be under- 
stood naturally. The horizon ra dius I is related to the 
cosmological constant A as I = ^/3/A. When I becomes 
large, A must decrease. Then the vacuum energy density 
A/87T also decreases. Thus, one finds intuitively that the 



true energy of the system never increases. The origin of 
the reversed sign for the energy lies on the negative bind- 
ing energy between matters and gravitational field of do 
Sitter space (35[. 

Returning to the calculation of the emission rate, we 
use dH = +duj and the replacement I — > l + uj. Then the 
imaginary part of the action becomes 



Im I = Im 



Jo r — {I + uj') + it 



drduj' '. 



(36) 



where we have used Feynman's prescription, lu' — > uj' 
Performing the integrations, eq. (|36|l results in 



-le. 



1 „ 

Im I = — 7T ( Iuj + — UJ 



(37) 



then the emission rate has the following form, 

r oc exp (2itIuj + ttuj 2 ) . (38) 

We find that the leading term gives the thermal Boltz- 
mann factor for the radiation, and that the second term 
represents the correction from the response of the back- 
ground geometry to the emission of a quantum. In fact, 
the leading term agrees with eq. fTS"]) . Hence, our inter- 
pretation about the energy in section[ll]is consistent with 
the result by the null geodesic method. That is, the lead- 
ing emission rate is written by means of [3c — 2irl and 
E = -uj as eq. (|27|) . 

If the entropy expression is used, equation (|3"5j) can be 
written in the more elegant form. Since the initial en- 
tropy of de Sitter space is Si = irl 2 and the final entropy 
is Sf = it(1+uj) 2 , one can write the emission rate, eq. 
as follows, 



T cx exp (S j - Si) = exp (ASbh) 



(39) 



This expression agrees with eq.([T]) which is the impor- 
tant result of ref.[5j. So, the emission rate is represented 
by the entropy difference between before and after the 
particle emission. 

Although we do not have the exact expression of the 
emission rate including the quantized gravity, quantum 
field theory tells us that the rate must be written as fl8l ] 



r cx \m 



(phase space factor) 



(40) 



where Mfi is the amplitude for the emission process. 
The phase space factor is obtained by summing over the 
final states and averaging over the initial states. Since 
the number of states is just the exponent of the entropy, 
the phase space factor is given by e Sf /e * = exp (ASbh)- 
This is nothing but eq. ([1} or (|39[) and agrees with the our 
result. Therefore our result shows that the emission rate 
deviates from the pure thermal emission and it is in full 
consistency with the unitary theory. 

Finally, we note that our consideration of the tunneling 
from the cosmological horizon have used the null geodesic 
method. In this method, since the positive frequency par- 
ticle follows its geodesic, the cosmological horizon must 
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expand or the final radius of the cosmological horizon 
must become larger than the initial radius. Suppose that 
the initial horizon is = Zj and the final horizon 77 = If. 
Then, since the horizon radius is related to the cosmolog- 
ical constant A as I = \/3/A, the relation Ti < rf leads 
us to the result 

A f <A l , (41) 

where Aj is the initial cosmological constant and A^- is the 
final one. This means that the value of the cosmological 
constant decreases by the particle emission from the cos- 
mological horizon. This will afford an active credence to 
our previous work on the thermodynamic consideration 
of the cosmological constant [30( ■ 

Our consideration can be applied to not only empty de 
Sitter solution but also black hole solutions in de Sitter 
space which have two event horizons. Some investigations 
similar to refs. [H, [l9[ have been done for these solutions 
[2(| HH [H, , in which they have considered the process 
that the cosmological horizon contracts. Thus, their re- 
sults may not describe the tunneling process by the null 
geodesic method on the correct way. For the cosmological 
horizon of black hole solutions in de Sitter space, simi- 
laly, one needs to consider the horizon expanding process 
in order to assure quantum tunneling. 

V. DISCUSSION 

In this paper, using the Hamilton- Jacobi and the null 
geodesic methods, we have investigated the emission rate 
from the cosmological horizon by quantum tunneling. We 
have shown that the leading order emission rate is given 
by ea. tfTij)) or (|27p and that the rate which includes the 
response of the background geometry is given by eq. (|38|) 
or P]). 

There are two important key points in our considera- 
tion. The first is the interpretation of the emission rate 
by the sign reversed energy, eq. (f2T|) . In de Sitter space, 
the energy of a black hole with mass m is evaluated as 
E = — m (34l ]. Along the similar line, we have proposed 
that the emitted particle should have energy E = —u>, 
where ui is the positive frequency of the particle. Then 
temperature of de Sitter space corresponds to the usual 
result (3c — 2ttI 321] . The other key point is the change of 
the radius of the cosmological horizon in the null geodesic 
method. When the particle's self-gravitation is taken into 
account, the horizon should expand in order to assure 
quantum tunneling. If the horizon contracts, the pair 
created particle can not tunnel out the horizon. The au- 
thors of refs. [H, have assumed that by a pair creation 
the background geometry changes from empty de Sitter 
space to Schwarzschild-de Sitter space. This, however, 
leads to the contraction of the cosmological horizon. Al- 
though they have considered that the positive frequency 
particle collapses to a black hole, the unobserved negative 
frequency particle does not contribute to the background 
geometry in their process. It contradicts with the original 



idea of quantum tunneling constructed in ref.Q. From 
this reason, we have reconsidered the tunneling process 
in this paper. 

Our study has an another important implication about 
the cosmological constant. When the radius of the cos- 
mological horizon expands, the value of the cosmological 
constant must decrease through the relation I = y^3/A. 
If A does not decrease, the cosmological horizon does not 
expand. Therefore, the pair created negative frequency 
particle should contribute to the background geometry 
as the horizon radius expands or the value of the cos- 
mological constant decreases. This will provide a strong 
confirmation to our previous work about thermodynamic 
consideration of the cosmological constant [3(| • 

The decrease of the cosmological constant has physical 
meaning as follows. First, the second law of thermody- 
namics holds for de Sitter space. Since the entropy of 
de Sitter space is given by S — 3n/A, S increases ac- 
cording to the decrease of A. Furthermore, the third law 
of thermodynamics also holds for de Sitter space. Tem- 
perature of de Sitter space is given by T# = l/2irl or 
Th = (l/27r)^/A/3. As quantum tunneling occurs, A 
continues to decrease but never be zero or negative. This 
is because the tunneling never happens if the final horizon 
does not exist. In order to exist the cosmological horizon, 
the cosmological constant A must be positive. Therefore, 
de Sitter space remains to be de Sitter space decreasing 
A gradually. According to the third law of thermody- 
namics, temperature of the system never reaches to zero 
by finite physical processes so that the cosmological con- 
stant never has the zero value. Even if A approaches zero, 
it never takes the exact zero value. Second, the Painleve 
coordinates (|3ip and planar coordinates 

ds 2 = -dt 2 + e 2t/ \dp 2 + p 2 dn 2 ) (42) 

are related each other by the coordinate transformation 
r = p e*' 1 so that the time coordinates are the same in 
both coordinates. Thus, the decay rates of the cosmolog- 
ical constant per unit time are also the same. Planar co- 
ordinates describe the exponential expansion of our uni- 
verse with a scale factor e l l l . Therefore, our universe will 
be de Sitter with decreasing A in future. From this point 
of view, our result does not seem to contradict with the 
observational value of the cosmological constant which is 
extremaly small. Third, if one believes the no boundary 
proposal for the quantum creation of universe by Hartle 
and Hawking [36[ , the most probable state is given by the 
nearly zero cosmological constant state [13, HH . The de- 
crease of the cosmological constant by quantum tunneling 
may imply that our universe gradually makes a transition 
towards the most probable state, that is, A — > state. 
Our result indicates that semiclassical consideration docs 
not contradict with these early works. 

The thermal decay of the cosmological constant has 
been studied recently in refs. [H, H(| • They have treated 
the decay process as quantum tunneling or thermal acti- 
vation. It is the natural next step to consider the relation 
between their decay process of the cosmological constant 



9 



and our tunneling approach for the cosmological horizon. radiation from the cosmological horizon, the underlying 
Although we do not still have the detailed mechanism physics is probably the same. We hope that this tanta- 
for the decay of the cosmological constant by Hawking lizing problem will be understood in the future. 
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